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ABSTRACT 

Pionium lifetime corrections are evaluated in the frameworks of constrained Bethe- 
Salpeter equation and chiral perturbation theory. Corrections of order 0(a) are calculated 
with respect to the conventional lowest-order formula, in which the strong interaction ampli- 
tude has been calculated to two-loop order with charged pion masses. The total correction 
is found to be of the order of (6.1 ± 2.0)%. 
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1. Introduction 

Our aim is to evaluate the sizable corrections to the nonrelativistic formula of 
pionium lifetimeB: 
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Here, Oq is the strong interaction (dimensionless) S-wave tttt scattering length in the 
isospin I channel and ^ + _(0) is the ground state wave function of the pionium at the origin 
(in x-space). The strong interaction scattering lengths a° and have been evaluated in 
^ . the literature in the framework of chiral perturbation theory (%PT) to two-loop order of 
the chiral effective lagrangianfl !' @. 

Formula ([!]) has been obtained by treating the strong interaction scattering am- 
plitude at threshold as a first-order perturbation with respect to the Coulomb potential 
between the charged pions. We shall evaluate in the following the corrections of order 
0(a) (a being the fine structure constant) to formula ([!]). Most of the details of the 
calculations can be found in Ref. i. 

2. Bound state equation 

To deal with the bound state problem, we use the constraint theory approach to the 
Bethe-Salpeter equation!, which consists in reducing the latter, by means of a constraint 
applied to the external relative energy, to a three-dimensional equation. The constraint 
that is chosen is the following: C(P,p) = (p\ ~Vt) ~ ( m i — m l) ~ 0> where p\ and P2 are 
the external particle momenta and m\ and m 2 their physical masses. One then establishes 
a three-dimensional eigenvalue equation that takes the form 

g-iy = -y*, (2) 
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where % l is the wave equation operator defined for two-spinless particle systems by 
g = —l/{p\ — m\ + ie)\ c = — \j{jp\ — m 2 , + ie)\ c , the index C denoting the use of 
the constraint. V is the potential, related to the renormalized off- mass shell scattering 
amplitude T by a Lippmann-Schwinger type equation: 

V = f + Vg f, f = -^=T, (3) 



2v^ 



c 



where s = P 2 = {pi + P2) 2 ■ The amplitude T contains the usual Feynman diagrams, 
where the external particles are submitted to the constraint C. The second term in the 
right-hand side of the first of Eqs. (|3|) generates an iteration series, the diagrams of which 
are called "constraint diagrams" , where the integrations relative to the factor g are three- 
dimensional, because of the presence of constraint C . The constraint diagrams cancel, 
in the s-channel, the singularities of the two-particle reducible diagrams of T. As long 
as perturbation theory is concerned, Eq. (|2]) is equivalent in content to the exact Bethe- 
Salpeter equation, with, however, a different arrangement of the perturbation series. The 
interest of Eq. (Q) is that it generates a three-dimensional systematic perturbation theory 
with covariant propagators!. 

We use for our calculations the chiral effective lagrangiani in the SU (2) x SU (2) 
case. We shall be working with the pseudoscalar densities defined as P a = iq^^q 
(a=l,2,3), where q are the quark fields and r a the Pauli matrices. The off-mass shell 
scattering amplitude is defined as the amputated connected four-point Green's function of 
the pseudoscalar densities, multiplied by the corresponding wave function renormalization 
factors 2B F Z~ 1 / 2 (p a ), where F is the pion decay constant F n in the chiral limit (F n = 
92.4 MeV) and Bq the quark condensate parameter!. With the above definitions, the 
amplitude T is renormalization group invariant. Furthermore, Eqs. (|^) imply that the 
constraint propagators g$ A b are kinematic operators with physical masses and are not 
concerned with the renormalization factors. 

For the analysis of the decay pfihe pionium into two neutral pions, it is convenient 
to use a coupled-channel formalism^' 13. Defining wave functions, potentials and propa- 
gators in matrix form with the indices {H — } for the charged pion sector and {00} for 
the neutral pion sector, Eq. (|2|) is transformed into two coupled equations. Eliminating 

the wave function \l/oo m terms of \E' H and neglecting unimportant unitarity correction 

factors, one ends up with a single wave equation for X I / + -: 



9o, 



Vcoul. + V+_ +_ — -^-,00,^0,00 ^00,- 



(4) 



represents the potential relative to the transition in the charged sector and from 
which the Coulomb potential Vcoui. = —m n +a/r has been separated. Voo,+- represents 
the potential relative to the transition from the charged to the neutral sector, etc.. In the 
constraint propagators <7o,+- an d <?o,oo h is the physical masses of the charged and neutral 
pions that appear, respectively. The nonrelativistic Coulomb potential will be considered 
as the zeroth-order potential and the remaining potentials in Eq. ([J) will be treated as 
perturbations. 

At the pionium energy, the operator ^ ,oo m Eq. (Q) lies in the scattering region of 
the 7T°7r sector and therefore contributes with imaginary values in x-space, thus inducing 
an imaginary part to the pionium energy. 
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3. Separation of quark-photon interaction terms 

In first-order of perturbation theory, it is the third term of the right-hand side of Eq. 
(|j) that contributes to the imaginary part of the energy and hence to the decay width. The 
simplest approximation at this stage consists in taking for the potential Voo,+- (= ^+-,oo) 
the strong interaction amplitude Tq^: at threshold (deviations from the pionium energy 
to the charged pion threshold yielding corrections of order a 2 ), calculated to two-loop 
order. Defining Pq = Pqr — iT/2, one recovers, at leading order in Am n , the nonrelativistic 
formula (|]). Corrections of order a to Tq arise from the presence of additional terms of 

electromagnetic origin in V+- t +- and Voo,h , as well as from second-order perturbation 

theory contributions of the potentials. It is natural at this stage to isolate from the 
rest the electromagnetic interactions coming from quark-photon interactions. The reason 
for this separation is that electromagnetic interactions with the explicit presence of the 
photon field lead, in bound state problems, or near threshold, to infra-red singularities 
which necessitate a detailed treatment. This is not the case of the quark-photon type 
interactions which operate only through mass shift effects and lead to regular expressions 
that can readily be continued from the bound state energy domain to the 7r + 7r~ threshold. 

The quark-photon interaction is represented in the chiral effective lagrangian at 
order e 2 p° i by the term e 2 C(QUQW) , where U is the chiral field, Q the quark charge 
matrix and C an unknown constant. Using for the field U the representation U = a + 
i7r.T/Fo, a = yjl — 7t 2 /Fq, where r are the Pauli matrices and 7r the pion fields, one 
finds: 

e 2 C(QUQU ] ) = -^tt+tt-. (5) 

This term induces for the charged pions the mass shift (Am 2 ) = (m 2 + — rrv^o) = 
2e 2 C / Fq , which is nonvanishing in the chiral limit and is numerically close to the physical 
pion mass shift! E. It has no effect on the scattering amplitude in lowest order, but acts 
essentially through insertions (together with counterterms of higher order) in the pion 
loop propagators, where the charged pion masses are replaced by their (almost) physical 
masses. 

We then split the scattering amplitude M. (= — iT) into two parts: 

M = M str - +qi + M em - , (6) 

where Ai ern " contains all electromagnetic interaction terms generated by explicit photon 
fields of the chiral effective lagrangian, while J\A str ' +q ' y contains the strong interaction 
terms together with the quark-photon interaction terms. 

Let us now neglect the amplitude Ai em \ It can be showni that the potential 
Voo,+-, calculated from Eq. ([|), where one takes also into account the contributions 
of the constraint diagrams, is equal, to order a, to the real part of the on-mass shell 
scattering amplitude Tqo^ 4 ! 97 taken at threshold: 

VS££P = Kef^(p 2 =p 2 2 = ml + ,p 2 3 = V \ = m 2 ^ s = 4ml + ,t = u = -Ami), (7) 

where Am 2 = m\+ — m\ Q . (An error in Ref. I concerning; the mass-shell condition of the 
neutral pions has been corrected.) 
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4. Second-order perturbation theory 

At this order, it is the interference of the last two potential terms of Eq. (d) that 
contributes to the imaginary part of the energy. 

The contribution of the discrete states of the pionium spectrum is found to be: 
(Ar)^ cr . = ( a /3)(2a° + a§)r . 

The contribution of the continuous spectrum is ultra-violet divergent, with a linear 
and a logarithmic divergences. In order to avoid these divergences, it is necessary to 
associate with the previous terms contributions coming from the constraint diagrams 
present in Eq. The constraint diagram with one pion loop has a linear divergence 
that cancels the previous one. The three constraint diagrams with one pion loop and 
one photon exchanged between the two internal pions have logarithmic divergences the 
sum of which cancels the previous logarithmic divergence. The sum of the contributions 
of the continuous spectrum and of the constraint diagrams is then finite and one finds: 

meant. + (^)const, = ^MK + 4)^ 

The total contribution of the discrete and continuous states, together with the 
constraint diagrams, is: 

(Ar) str = 1.5a(2o§ + a 2 )T = 0.004r . (8) 



5. Vacuum polarization 

Vacuum polarization in positronium and hydrogen-like atoms provides corrections 
of order a 3 . In the pionium problem, vacuum polarization, because of the presence of 
electron loops and of the large ratio of the reduced mass of the two-pion system to the 
electron mass, yields corrections of order a. This point was emphasized by Labelle and 
Buckley0. Using the nonrelativistic QED effective field theory formalism^!, they have 
found: (Ar) mc pol = 0.43oT . 

In the present formalism, vacuum polarization contributes through the potential 
V +-.+- with a local termlil. It contributes to the decay width through the interference 
term between V + _ )+ _ and the last term of Eq. (^) at the second order of perturbation 
theory. The contribution of the discrete states is found to be: (AT) discr = 0.03ar , 
while that of the continuous states, which is finite, is (AT) cont = 0.38ar . The total 
contribution is: 

(ATL vol. = 0-41aT = 0.003r , (9) 
in agreement with the previous estimate. 



6. Electromagnetic radiative corrections 

These corrections arise from pion-photon interaction terms contained in _M em " [Eq. 
(§)]. The chiral effective lagrangian for the SU(3)xSU(3) case including electromagnetism 
to order e 2 p 2 was presented by Urechi. It was presented for the SU (2) x SU (2) case by 
Meissner, Miiller and Steininge and by Knecht and UrechEl. We use for our subsequent 
calculations the SU (2) x SU (2) lagrangian with the notations of Ref. 

The pion-photon radiative corrections to order e 2 p 2 are represented by the pion 
self-energy corrections and the vertex correction of the lowest order scattering amplitude 
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of the process tt + 7t~ — > 7r°7r°. The sum of these contributions to the off-mass shell scat- 
tering amplitude M. m+ _ at the pionium energy has a spurious 0(a°) term. The latter 
is cancelled, however, by the contribution of the corresponding constraint diagram of Eq. 
(i) ("with one photon exchanged between the charged pions). The corresponding total 
contribution is: 



KAem. 
- /Vt 00,+ 



» jem.iconstr.) 
+ - /Vl 00,+- 
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^0 



,(10) 



where the k r, s are the renormalized low energy constants appearing in the counterterm 
lagrangian and /i is the renormalization mass, chosen in the following for the numerical 
estimates equal to the p- meson mass (770 MeV). We have also incorporated in the mass 
term of the tree level amplitude the mass shift of the neutral pion. For the numerical 
estimate of the corrections, we have adopted the values of the k r, s obtained by Baur and 
UrechU (in the Feynman gauge) using a resonance model for the saturation of sum rules 
and assigned to them conservative 100% uncertainties. The numerical values that we use 
(after conversion to the SU(2) x SU(2) case) are, in units of 10~ 3 : k T 2 + 10/c[ = 122.2, 
/C3 = 6.4, kj + 2fcg = 0.8. One finds for the modification of the decay width: 



(Ar)^ 7 = ( - 0.0015 ± 0.0075)r c 



where the uncertainty comes from the fc r, s. 



7. Electromagnetic mass shift corrections 

These corrections are those contained in the amplitude 7£eA^ocvf- 7 [Eqs. ©-©]• 
They are induced by the term (j5|) and by its counterterms of the effective lagrangian 
and the effect of which is the generation of mass shifts in the scattering amplitude. The 
corrections are evaluated with respect to the strong interaction scattering amplitude, in 
which all pions have the same mass, chosen to be the charged pion massa H @. 

In the strong interaction case, the pion mass is essentially determined in terms of 
the quark mass parameter m a and the quark condensate parameter Bq. In lowest order, 
one has m\ = 2mBo. To one-loop order, one has: 



(mir- = 2mB + + J-jlnfe) = 2mB, (12) 



where 1% is one of the low-energy constants of the effective lagrangian. 

According to the result @, the neutral pion mass remains unchanged when the 
0(e 2 p°) term is introduced in the strong interaction lagrangian. Therefore, to order p A 
of the strong interaction lagrangian and to order e 2 p° of the electromagnetic interaction 
the mass appearing in the left-hand side of Eq. ( |l~2|) can be identified with the neutral 
pion mass. Higher-order effects in the masses play only a corrective role and modify 
slightly this equality, which can be represented in the approximate form as m 2 = 2mB + 
0(e 2 p 2 ) + 0(p 6 ) ~ 2mB. 
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Let us now rewrite the amplitude IZeAi^'+T 1 by specifying in detail its kinematic 
and parametric conditions, as emerging from Eqs. (|7|) and (0): 

KaM$+r = KeM^-i* = Aml + ,pl = p\ = = p\ = 2mS ~ m\ a ). (13) 

On the other hand, the strong interaction amplitude is defined in the following way: 

KeM^ ;+ _ = KeM^; + _(s = 4m^,p? = p\ = ml +:V \ =v\ = ml + ,2mB = m 2 + ). (14) 

We have to calculate the difference between the two previous amplitudes. The latter is 
the result of three effects: (i) The dynamical effect of the quark-photon interaction term 
([J), which acts through insertions in internal pion loop propagators, (ii) Shift of the 
neutral pion momenta squared, p\ and p 2 , from their value m 2 + of the strong interaction 
case to their final value m 2 ^. (iii) Shift of the mass parameter 2mB from (strong 
interaction case) to m 2 ^. As long as one considers linear effects in the pion mass shift, 
these three effects can be calculated separatelyi. We shall stick to this approximation in 
the following. One then obtains: 



[ATZeM 



Ami mlo Ami f 22 ,fnlo N 23 



F 2 16tt 2 F 4 



jj 2 6 

(^2 2^) ~~\~ \k"J 



(15) 



The T's and k r, s are the renormalized low energy constants of the strong interaction and 
electromagnetic interaction lagrangians, respectively. The logarithms that appear in the 
strong interaction amplitude, associated with the T's, are calculated with the charged pion 
mass. The masses that appear in the corrective terms are taken equal to the neutral pion 
mass. The Am 2 factor of the lowest-order term represents the physical pion mass squared 
shift, for the 0(e 2 p 2 ) terms of the neutral pion mass, coming from pion-photon and quark- 
photon interactions, have been included in the tree level amplitude. We neglect, in our 
numerical evaluations, isospin breaking effects coming from the quark masses, which are 
smali3 m We use, as in Sec. 0, the numerical values of the /c r 's obtained in Ref. (in 
the Feynman gauge). These are, in units of 10~ 3 : k r 2 = 5.4, k\ = 6.2, k r 7 — 2k r 8 = —2.3, 
kj + k r xl = —0.3. The constants V are obtained from Refs. ; they are, in units of 10 : 
l\ = -5.4 ± 3.9, 11 = 0.8 ± 3.8, l\ = 5.6 ± 5.7. We also take F = 88 MeV I. 
The correction for the decay width is: 



(^r)m. shift 



0.049 ±0.004 ±0.001, (16) 



where the first uncertainty comes from the T's and the second one from the k r, s. The 
contribution to the above number of the tree level correction, represented by the first 
term of the right-hand side of Eq. (|U|) and coming from the effect of type (iii), is 0.037. 
Therefore, the 0(e 2 p 2 ) terms contribute with an amount of 30% with respect to the tree 
level effect. 
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Knecht and Urech!i3 have provided, from a direct evaluation of the generating 
functional to one loop, the expression of the on-mass shell 7t + tt~ — > 7r°7r° scattering 
amplitude in the presence of electromagnetism. A finite result at threshold is obtained 
after subtraction of the infra-red singularity of the photon. The sum of our evaluations 
( PD and ( p~5D agrees with the result presented by the above authors, to linear effects in 
Am n . 

8. 0(e 2 p 4 ) effects^ 

0(e 2 p A ) effects come from diagrams having a pion loop with one photon exchanged 
or having two pion loops with a mass shift insertion. Most of these effects are expected 
to be small. However, enhancements may appear due to the presence of infra-red sin- 
gularities. The latter generally arise from kinematic regions where the loop momentum 
k is small. Two cases are distinguished. In the first case, the component k is of the 
order of k 2 /m n ; in the second case, ko is of the order of |k|. The singularities of the first 
case, which are the strongest ones, are generally cancelled by those of the corresponding 
constraint diagrams. The singularities of the second case are generally cancelled either 
by those of the crossed diagrams (if present) or by those of the self-energy diagrams. A 
typical example of such cancellations was met in the vertex correction evaluation [Sec. 
EX where the sum of the contributions of the vertex diagram and of the corresponding 
constraint and self-energy diagrams was found to be free of singularities. 

By considering groups of diagrams of the above kinds, one arrives at the conclusion 
that the sum of all 0(e 2 p A ) diagrams is free of infra-red singularities (at the pionium 
energy). However, in our treatment of the second-order perturbation theory problem [Sec. 
[D, contributions of the constraint diagrams of the diagram with a photon exchanged 
between the two loop-pions were already taken into account to cancel the ultra-violet 
divergence of the continuous spectrum of the intermediate states. Therefore, the infra- 
red contribution of the corresponding four-dimensional diagram should also be isolated. 
This diagram has an infra-red singularity of the first kind mentioned above and has been 
evaluated in the literature Its contribution to the pionium decay width, leaving aside 
its ultra-violet divergent part, is: 

(Ar) G(eV) = ~{2a° + a 2 ) (2 In a + 3 In 2 + 21C(3)/(27r 2 ))r = 0.006r . (17) 

The remaining groups of diagrams, being infra-red regular, will have orders of 
magnitude fixed by x?T. Let us consider the 0(e 2 p 2 ) radiative correction of Sec. |6]. 
The corresponding relative correction to the decay width, without the counterterms k r , 
is of order a [Eq. (0)]. Pion loops have typical relative order of magnitude of 5% 
(compared to 1). The number of the 0(e 2 p A ) groups of diagram being of the order of 25 
and assuming their signs are distributed at random, one arrives at an estimate of 0.2% for 
the corresponding relative contribution to the decay with. A similar contribution, 0.2%, 
should also be expected from the low energy constants of the counterterm lagrangian. 
Mass shift corrections due to quark-photon interactions have produced a 30% effect at 
the 0(e 2 p 2 ) level relative to the tree level, with a 1% global effect [Eq. (|TB|) and comment 
following it]. Assuming that the corresponding 0(e 2 p 4 ) effect (two-loop diagrams) is of 

"'The content of this section has been added after the Workshop. I thank E. A. Kuraev, P. Labelle and A. G. 
Rusetsky for drawing my attention to this question and for discussions. 
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relative order of 40% with respect to the 0(e 2 p 2 ) effect (taking into account the increase 
in the number of diagrams), one arrives at an estimate of 0.4% for its contribution. The 
sum of the previous contributions is then 0.8%. The total contribution of the 0(e 2 p 4 ) 
diagrams is thus: 

(Ar) 0(e 2 p 4) = (0.006 ± 0.008)r . (18) 

9. Summary and concluding remarks 

The corrections to the nonrelativistic formula of the pionium decay width can be 
represented in the following form: 




where Ar is composed of the following contributions: 

Ar = (AT)^. + (AT) mc . poL + (AT) W1 + (AT) m shlft + (Ar) 0(eV) . (20) 
One finds for Ar/IV 

Ar , . 

— = 0.061 ±0.004 ±0.008 ±0.008, (21) 
To 

the first uncertainty coming from the low energy constants /[, the second one from the 
fej's and the third one from the 0(e 2 p 4 ) diagrams. T [Eq. (|l|)] is calculated with the 
strong interaction scattering lengths evaluated up to two-loop order (a° — a 2 , = 0.258 
!). One finds r = 3.19 x 10 -15 s. With the corrective terms, the lifetime becomes: 
t = (3.02 ±0.07) x 10- 15 s. 

The different contributions are summarized in Table 1. 



Effect 


Ar/r (%) 


strong 
vac. pol. 
7T7 rad. corr. 
mass shift 
O(eV) 


0.4 
0.3 
-0.1 ±0.7 
4.9 ±0.4 ±0.1 
0.6 ±0.8 


total 


6.1 ±2.0 



Table 1: The various corrections to the decay width. 

Evaluations similar to those of the present work were done by Ivanov et a/.IH, 
solving the Bethe-Salpeter equation in the Coulomb gauge and using, for the mass shift 
corrections, the result provided in Ref. El. The results found by Ivanov et al. and 
ours are qualitatively and quantitavely compatible among themselves. In particular, the 
sum of their evaluations corresponding to strong and one Coulomb exchange corrections 
(—0.22 + 1.55) should be compared with the sum of our evaluations of strong and 0(e 2 p 4 ) 
effects (0.4 ±0.6). The latter effect, not evaluated in Ref. i, has reintroduced an infra-red 
logarithm present in Ref. 111. The sum of the quantities corresponding to mass shift and 
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electromagnetic radiative corrections of Ref. li^l (2.99 + 1.73) should be compared with 
the sum of our evaluations of mass shift and pion-photon radiative corrections (4.9 — 0.1). 
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